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Abstract 

The transform considered in the paper averages a function sup- 
ported in a bah in over ah spheres centered at the boundary of 
the ball. This Radon type transform arises in several contemporary 
applications, e.g. in thermoacoustic tomography and sonar and radar 
imaging. Range descriptions for such transforms are important in all 
these areas, for instance when dealing with incomplete data, error 
correction, and other issues. Four different types of complete range 
descriptions are provided, some of which also suggest inversion pro- 
cedures. Necessity of three of these (appropriately formulated) condi- 
tions holds also in general domains, while the complete discussion of 
the case of general domains would require another publication. 
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1 Introduction 

The spherical mean Radon transform, which integrates a function over all 
spheres centered at points of a given set, has been studied for quite a while 
in relation to PDE problems (e.g., ^31^ EH]). However, there has been a 
recent surge in its studies, due to demands of manifold applications. These 
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include, among others, the recently developed thermoacoustic and photoa- 
coustic tomography (e.g., [131 E21 IS21 ES]; |^- |lU3j ). as well as radar and 
sonar imaging, approximation theory, mathematical physics, and other areas 
[ilZlliaillliHlEliniliniinSllZil^En]. For instance, in thermoacoustic 
(and photoacoustic) tomography, the spherical mean data of an unknown 
function (the radiofrequency energy absorption coefficient) is measured by 
transducers, and the imaging problem is to invert that transform (e.g., jH2], 
P^- |lU2j ). These applications also brought about mathematical problems 
that had not been studied before. Many issues of uniqueness and stabil- 
ity of reconstruction, inversion formulas, incomplete data problems, etc., are 
still unresolved, in spite of a substantial body of research available (e.g., 

isiiziii3iiinii2siiz3Enii2niE2iEniii3iMiinziisa 

EOllHaiHSllHllHHllHHllHniinni, inni-CnSl)- in this text we address the problem 
which has been recently considered for the first time ^3 [MJ |Hn] (see also 
related discussions in [73J [73] ) , namely the range conditions for the spherical 
mean transform. In fact, as we will mention below, in an implicit form, a 
part of range conditions was already present in |HTl lU^ and later in [7j. 

For someone coming from PDEs and mathematical physics, the range de- 
scription question might seem somewhat unusual. However, it is well known 
in the areas of integral geometry and tomography that range descriptions 
are of crucial theoretical and practical importance |2Z1 IS31 ISSl tHHl 1331 HOI IZ21 ■ 
The ranges of Radon type transforms usually have infinite co-dimension (e.g., 
in spaces of smooth functions, or in appropriate Sobolev scales), and thus 
infinitely many range conditions appear. One might wonder, what is the im- 
portance of knowing the range conditions. The answer is that, besides their 
analytic usefulness for understanding the transform, they have been used for 
a variety of purposes in tomography (as well as in radiation therapy planning 
[2211231 Eni): completing incomplete data, correcting measurement errors and 
hardware imperfections, recovering unknown parameters of the medium, etc. 
[36, £6, IHHl ISni HOI [Za [7H1 [HIl ins ESI- Thus, as soon as the spherical mean 
transform started attracting a lot of attention, researchers started looking 
for range descriptions. Some range conditions (albeit they were not called 
this way) were already present in [^H [n21 and in [7] (see also ) , where the 
sequence of polynomials was considered arising as moments of the spherical 
Radon data. In an explicit form, these conditions were formulated recently 
in [Sn]. They, however, as it was discovered in 15 , do not describe the range 
completely. A complete set of conditions was found in the two-dimensional 
case in the recent paper ^Hj and for odd dimensions (albeit, for somewhat 



3 



different transforms) in [31]. In all these papers, the centers of spheres of 
integration (i.e., the location of tomographic transducers) were assumed to 
belong to a sphere. 

In this paper, we obtain range descriptions in arbitrary dimension for 
the case of centers on a sphere. Moreover, we obtain several different range 
descriptions that shed new light on the meaning of the range conditions 
(in particular, onto the appearance of two seemingly different subsets of 
conditions) . 

In Section 121 we introduce main notations and preliminary facts that will 
be needed further on. Section El contains the formulation of the main results. 
Theorem ^1 provides three different types of range descriptions. Theorem 
^2 establishes that in odd dimensions moment conditions of Theorem ITUl can 
be dropped (a situation analogous to the one in |31i)- On the other hand, it 
is shown in Theorem ^1 that a strengthened version of moment conditions 
alone describes the range. It is shown in Theorem ^1 that the results also 
hold in Sobolev scale, rather than in category. The next four sections 
are devoted to the proofs of these theorems. It is noticed in Section |H1 that 
necessity of most of the range conditions is in fact proven for general domains, 
not just for a ball. This is described in Theorem|2H Section El contains proofs 
or alternative proofs of some technical lemmas. The alternative proofs are 
provided, since the authors believe that they might shed extra light onto the 
problem. The final two sections contain additional remarks and discussions, 
and acknowledgments. 

2 Main notions and preliminary information 

In this section we introduce main notions and notations that will be used 
throughout the paper. We also remind the reader of some known facts that 
we will need to use. 

We will be dealing with domains in M."^. The closure of a domain Q is 
denoted by Q and its boundary is dQ. We denote the n— dimensional unit ball 
-B = {x G M"! |x| < 1} and the unit sphere is S" = dB = {x G ]R"| \x\ = 1}. 
The area of S" is and the area measure on S will be dS (this notation will also 
be used for the surface measure on the boundaries of other domains). The 
notation C^{B) stands for the class of smooth functions with the compact 
support in the closed unit ball. For partial derivatives, the notations dtf, 
and ft will be used. 
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2.1 Spherical means 



The main object of study in this paper is the spherical mean transform R 
(with centers on S* = dB) that takes any function / G C^{B) to 

Rfip,t):=^ J f(p + ty)dSiy),peS. (1) 

One might wonder why we require the support of / to belong to B. It will 
be explained in Section El that there is not much hope for explicit range 
descriptions, if one allows the support of the function to spill outside the 
surface S of the centers. 

We will also consider the cylinder C = i? x [0, 2] and its lateral boundary 
S X [0,2]. 

A t 




Figure 1: The geometry of domains 



Notice that the sphere S enters here in two different ways: as the set of 
centers of the spheres of integration {p G S), and as a parametrization of 
the spheres of integration {Jy^g ...dS{y)). The reader should not confuse the 
two, since sometimes in the text we will change the set of centers to a more 
general surface S, while keeping the integration surfaces spherical. 

The results can be easily re-scaled to the case when the set of centers 
is a sphere of an arbitrary radius p. We avoid doing this here, in order to 
simplify the expressions. 
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2.2 Darboux equation 

Allowing in the centers x of the sphere of integration be arbitrary, one 
arrives to a function 

G{x,t) = - [ f{x + ty)dS{y),xeW. (2) 

y&S 

It is well known [121 1211 11H| that the function G{x, t) defined by satisfies 
the Darboux (Euler-Poisson-Darboux) equation 

as well as the initial conditions 

dG 

G{x, 0) = fix), — (x, 0) = 0, X G M". (4) 

Moreover, any such solution of (jH)) - © in M" x is representable as the 
spherical mean Q of f{x) (Asgeirsson's Theorem, see [T7 t 1241 1^ )- 

An important remark about the initial conditions Q is that they mean 
that the solution G can be extended to all values of time as an even solution 
on X 1^ Ch. VI. 13]. 

One notices that the restriction of G to S* x ]R+ coincides with Rf: 

g = GUxR^ = Rf (5) 

Another observation concerning the mean G{x,t) of a function f{x) sup- 
ported in B is that it vanishes for x E B,t > 2. Indeed, the value G{x,t) is 
the average of / over the sphere centered at x and of radius t, while such a 
sphere for x E B,t > 2 does not intersect the support of /, which is contained 
in B. So, G satisfies the terminal conditions 

dG 

G(x,2) = 0,— (x,2)=0,xGi?. (6) 

We also need to mention some other known properties of the Darboux 
equation, which we will need to utilize further on in the text. 

The Darboux equation has a useful connection with the wave equation 
(e.g., Ch. VI. 13]). This relation comes from existence of transformations 
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that intertwine the second derivative operator with the Bessel operator 
Bp = ^ + ^^r^^- ^ general approach to constructions of such transformation 
operators (not just for the Bessel case) can be found, for instance in jnO] (see 
also 51 ). Among those, the most commonly used ones are the Poisson 
transform (also called Delsarte or Riemann-Liouville transform [211 123 
inSl HZ] ) and Sonine transform ^ . The formula we will need for the Poisson 
transform is 

t 

iV.U) it) = c,t-^^ [ U{y){e - yy-'''dy. (7) 



Here Cp is a non-zero constant, whose specific value is of no relevance to us. 
We will use this transform for even functions f/, in which case it can be 
rewritten as 

1 

{VpU) (t) =CpJ Uifxt) (1 - fxy-'/'dfi. (8) 

The properties of this transform are well known (e.g., |^ EHl IHOl EZI)- In 

particular, it is known to preserve evenness of functions. We will also need 
the inversion formula in the case when 2p is an odd integer: 

{V-'G) it) = const t j [t'^Git)) . (9) 

We will not use the Sonine transform here. We, however, will be interested 
in the intertwining operator sometimes called Weyl transform^ p7] : 

Uit) = {WpG){t) := Jl^^l^^^^^ I G{s){s'-ty-yhds. (10) 

t 

which we will use for the specific value p = {n — 2)/2. The inverse transform 
W^^ for this case is given by 



oo 



Git) = (w-'u) (t) = { ^ 



const / (s^ — t^) ^/^(^^^)"/^f/(s)s(is, for n even 

)dd. 

(11) 



const^T^ ^ U{t)^ for n odd. 



-'^AU such transforms that we use are partial cases of transforms of Abel type. 
^This transform is dual to a Poisson transform |51l I97j . 
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One has (e.g., [HIl EHl EZl ) the following intertwining identities: 



These identities show that G{x, t) (for t > 0) is a solution of Darboux equa- 
tion if and only if its Weyl transform (with respect to t) U (x, t) = WpG{x, t) 
solves the wave equation. 

It will be important for us that the transform W (unlike the Poisson 
transform) involves integration from t to oo. Hence, when applied to a func- 
tion on that vanishes for t > a, it preserves this property. It is clear from 
the inversion formulas that behaves the same way. 

An important relation between Fourier, Fourier- Bessel, and Weyl trans- 
forms will be presented in the next subsection. 



= (|i + ^1) 



2.3 Bessel functions and eigenf unctions of the Bessel 
operator and Laplacian 

It will be convenient for us to use a version of Bessel functions that is 
sometimes called normalized, and sometimes s;?/ienca/ Bessel functions (e.g., 

isnna ?]): 

. !!£(£±im. (13) 

For p = according to Poisson's representations of Bessel functions 
IU7] . this is just the spherical average of a plane wave in M". 
We will also use standard expansions 

A- 



(14) 



JpW 2P ^( -'-) 22'=fc!r(p+A:+l)' 
fc=0 

oo 
k=0 

where Ck are non-zero constants. 

The function z{t) = jp{Xt) satisfies the Bessel equation 

d^z 2p+ldz 2 /irN 
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and initial conditions 

z(o) = i,4(o) = o. 

The standard Fourier- Hankel (also called Hankel, or Fourier-Bessel) trans- 
form g{t) ^ J^p{g){sX) and its inverse can be nicely written in terms of jp-. 

oo 

00° (16) 
9{t) = ' J J'p{g){tX)UXt)X'P+'dX. 


We will use notation JF for the standard one-dimensional Fourier trans- 
form. 

The following relation between the Fourier, Fourier-Hankel, and Weyl 
transform (e.g., [HZI P- 124]) helps to understand some parts of the further 
calculations: 

J^p = CpTY^p, (17) 

where Cp is a non-zero constant, explicit value of which is of no relevance to 
this study. This relation shows that the Weyl transform is the ratio of the 
Fourier and Fourier-Bessel transforms. 

We will need to use the following known Paley- Wiener theorem EHl 
IMJinZI for Fourier-Bessel transform^. 

Lemma 1 Let-p = (n— 2)/2 and H be a Hilbert space. An H -valued function 
^(A) on M can be represented as the transform 17^) of an even function 
g G C^{W,H) supported on [—a, a], if and only if the following conditions 
are satisfied: 

1. $(A) is even. 

2. $(A) extends to an entire function in C with Paley-Wiener estimates 

||$(A)|| < Cn{1 + |A|)-^e'^l^"^ ^1 (18) 

for any natural N. 

■^Albeit the statement of the Lemma is formulated in the works cited only for the case 
H = C, the proofs allow one to consider without any change valued functions as well. 
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Although this result is well known, for reader's convenience, we provide its 
proof in Section El 

We will need some special solutions of Darboux equation in the cylinder 
B xM.. Let — be in the spectrum of the Dirichlet Laplacian in the ball B, 
and ip\{x) be the corresponding eigenfunction, i.e. 

Ai^xix) = -X'^tpxix) in B ^^^^ 
ipx{x) = on S*. 

Equations (fT^ and (fT^ imply that the function 

MA(x,t) = V^A(a;)jn/2-i(At). (20) 

satisfies Darboux equation Q. 

Finally, we need descriptions of (generalized) eigenfunctions 

-A^ = X'^tp. 

of the Laplace operator in M"and in B in terms of their spherical harmonics 
expansions. 

Let Y[^{6), 6 G 5*""^, m = 0, 1, 1 < / < d{m), be an orthonormal basis 
of spherical harmonics, where m is the degree of the harmonic and 

T{m + n — 2) 

d[m) = [2m + n — 2)- 



r(n- i)r(m + !)■ 



It is known (and can be easily shown by separation of variables) that any 
function 

<pmA^) = (Ar)i-"/V„/2-i+™(Ar)l^-(^) , 

= (Ar)'"j„/2-i+m(Ar)y;'"(^), where x = rO ^ ' 

is a generalized eigenfunction of the Laplace operator^ in M" with the eigen- 
value — A^. In fact, one can show that any generalized eigenfunction of the 
Laplace operator in M" has the following expansion into spherical harmonics: 

u{re) = 5^Q,„(Ar)"^j„/2-i+^(Ar)F™(^), (22) 

Lm 



^ As it is customary, we use the term "generalized eigenfunction" for any solution of the 
equation Au — ~\^u in M". 
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where r = 1x1,6' = - — r. One can describe precisely the conditions on the 

\x\ 

coefficients q^^, necessary and sufficient for to provide all generalized 
eigenfunction, as well as generalized eigenfunctions with some prescribed 
growth condition at infinity pi |2] . 

If one chooses only the values of A 7^ that are zeros of Jm+n/2-iW , one 
arrives to the eigenfunctions of the Dirichlet Laplacian in the unit ball B. In 
particular, functions (PT|) for all m,l and A 7^ such that Jm+„/2-i(A) = 0, 
form a complete set of eigenfunctions. 

With all these preparations in place, we can now set out to formulate and 
prove the results of this article. 



2.4 Some preliminary results 

We start considering the spherical mean transform R introduced in (^J with 
centers on the unit sphere S. The following moment conditions for this 
case were present in IHH El], as well as in [7j (see also [2QJ), and explicitly 
formulated as range conditions in 



Lemma 2 Let g{x,t) = Rf{x,t) for f G C^{B). Then, for any non- 
negative integer k, the function Mk on S defined as 



2 



Mkix) := j t^''+''-^g{x,t)dt,x E S (23) 







has an extension to M" as a polynomial qk{x) of degree at most 2k. 
Proof One readily observes that, for |x| = 1, 

Mk{x) = j \x-p\''f{p)dp. (24) 

Applying ()24j) to arbitrary x G M" (not necessarily on the unit sphere), one 
clearly gets a polynomial Qkix) of x of degree at most (not necessarily equal 
to) 2k. □ 

Remark 3 In fact, noticing that 

Mk{x):= [ \x-p\^^f{p)dp= f{l-2p-x+\p\^ff{p)dp, 
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we conclude that Mk{x) has an extension to as a polynomial of degree at 
most k. 

The possibility of reducing to degree k comes from the fact that we use 
centers of the spheres of integration that belong to the (unit) sphere. It is 
clear that when the centers run over a non-spherical surface, this reduction is 
not possible anymore, and one cannot guarantee degree less than 2k. How- 
ever, in this case an extra condition can be found that straightens up the 
situation. We thus provide here an alternative reformulation of the moment 
conditions, which will be handy in more general considerations. 

Lemma 4 Let D G be a smooth bounded domain and R denote the 
spherical mean transform with centers on S = dD. 

Let g{x, t) = Rf{x, t) for f G C^{D). Then, for any non-negative integer 
k, the function Mk on S defined as in \2'^) has an extension Qk{x) to M" as a 
polynomial of degree at most 2k, satisfying the following additional condition: 

^Qk = CkQk-i for any k = 1, (25) 
where Ck = 2k{2k + n — 2). 

Proof Let us notice that = jxp'^ * f{x)\s- We can now define the 
polynomials Qk{x) = {x]"^^ * f{x). Then the relation follows from the easily 
verifiable identity A|a;p'^ = Ck\x\'^^''~^\ □ 

In fact, in the case when D is a ball B, the condition is not needed. 

Lemma 5 In the case when D = B, the condition \25]) of Lemma ^ on a 
function g{x,t) on S x [0,2] can be dropped, and thus conditions of Lemmas 
\^and^are equivalent. 

Proof To prove the lemma, we need the following well known fact (e.g., |49j). 
which we prove here for the sake of completeness. 

Proposition 6 The solution of the boundary value problem 

Au = V, \x\ < 1, 
u = 0, \x\ = 1, 

where v is a polynomial, is a polynomial of degree degu = degv + 2. 
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Proof Let us prove first that there exists a polynomial solution u of Poisson 
equation 

Au = V 

in the unit ball, such that degu = degv + 2. Clearly, it suffices to do this 
for each homogeneous term of f , so we can assume the polynomial v to be 
homogeneous. 

Let us represent v in the form: 

iy=0 

where each hi, is either zero, or a homogeneous harmonic polynomial of de- 
gree deghu = degv — 2z/, and brackets [...] denote the integer part. This 
representation is well known to be always possible (e.g., [13 EE])- A solution 
u can be found in the similar form (where we denote for brevity k = [^^^^j^]) 

k 

u{x) = \x\'^'^hu{x) , deg hy = degv + 2 — 2v. 

Here again, each hjy is either zero, or a homogeneous harmonic polynomial 
of degree deg h^, = deg v + 2 — 2z/. Then direct calculation shows 

k 

Au{x) = ^ (^Cu + 2i/deghjj \xf'^~'^hi,{x), 

iy=0 

with the coefficients c^, defined in Lemma 3. Thus, the needed polynomial 
solution u that we are looking for can be obtained by choosing 

hy = [cy + 2v{degv + 2 - 2v)Y^K^i. 

To finish the proof of Proposition, introduce U = u — u. Then one obtains 

AU = 0, |x| < 1, 

U = u, \x\ = 1, 

for the newly defined function U. The boundary value u is the polynomial 
of degree degtt = degf + 2. Its harmonic extension U from the unit sphere 
is obtained from the above decomposition of u by replacing by 1: 

k 

U{x)=^K{x). 

iy=0 

13 



Since deg h,^ < deg f + 2, we have deg U < deg v + 2. 

Thus, the solution u = u — U is a polynomial of degree at most degv + 2, 
which proves the Proposition. □ 

Remark 7 Polynomial solvability of the Poisson problem with polynomial 
data is rather unique and essentially holds only for balls (e.g., 1211 \^). 
Hence, in a general domain, the conditions of Lemma ^ are stronger than 
the ones of Lemma{^ 

Proof of Lemma We want to prove that among all polynomial 
extensions qk, degqk < 2k, of the functions Mk defined in (17), there is a 
sequence of extensions Qk satisfying the additional recurrence relation (j2S)). 
Let qk be some extensions. Any other sequence Qk of extensions can be 
represented as 

Qk = Qk + Uk, 

where = on the unit sphere. The additional requirement ((23) yields the 
relation 

Auk = CkQk-i - Aqk. 

The existence of polynomial solutions Uk, deg Uk < 2k follows now by induc- 
tive application of Proposition IHl Then the modified sequence of polynomials 
Qk = Qk + Uk satisfies all the requirements of the lemma. □ 

Remark 8 The sequence Qk of polynomial extensions of the functions Mk 
satisfying the chain relation (19) is unique. Indeed, if Q'k is another sequence 
of such extensions, then the polynomials Rk = Qk ~ Q'k vanish on the unit 
sphere and still possess \25]) . If not all Rk are identically zero, then, due 
to \25\) . the first nonzero polynomial Rk^ is harmonic and vanishes on the 
unit sphere. Thus, it must be zero, due to the maximum principle. This 
contradiction shows that Rk = and hence Qk = Q'k for all k. 

In the future, we will need the moment condition on the unit sphere in a 
different form |15j : 

Lemma 9 The function g{x,t),x G 5", t G [0,2] satisfies / f^) if and only if, 
for any spherical harmonic ¥"^{6) of degree m, the function 

2 

gmW = j j g{e,t)j^/2-i+n.{\t)Y"\e)dS{e)e-^dt 

ees 

has at X = a zero of order at least m. 



14 



Proof The moment conditions require that 

t'^+''-^g{e,t)dt 







is extendable to a polynomial of degree at most k. Let us expand g{6, t) into 
an orthonormal basis Y^^ of spherical harmonics on S (where m is the degree 
of the harmonic): 

g{e,t) = Y,9iAt)yno). 

l,m 

Due to smoothness and compactness of support of g, it is legitimate to inte- 
grate term-wise in computing the momenta to obtain 

oo oo 

j f'+'''-'g{9, t)dt = J2 Yrid) J t'''^''-'9i,Mdt. (26) 



A spherical harmonic of degree m can be extended to a polynomial of degree 
d if and only if d > m and — m is even. Thus, the moment conditions 
require that the coefficients t^^+''-^gi,m{t)dt in ^ with m > 2k must 
vanish. 

Let us turn to the function c/mW- Using the expansion (fT^ for Bessel 
functions, one arrives to the formula 

CO p p 

gUX) = J2 CkX'" J J Y"'{e)g{e, t)ds{e)dt. (27) 

Now one sees that the moment conditions are equivalent to the requirement 
that all terms in this series with 2k < m vanish. Therefore, the series begins 
with the power at least A"^. □ 



3 Statements of the main results 

As it follows from |^ , the moment range conditions of the preceding 
lemmas are insufficient. Necessary and sufficient conditions of different kinds 
were provided in ^3] in dimension two and in |31j in odd dimensions (albeit 
for somewhat different transforms). We formulate below our main result 
that resolves the problem of range description in any dimension, as well as 
provides several alternative ways to describe the range. 
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Theorem 10 The following four statements are equivalent: 

1. The function g G C^{S x [0,2]) is representable as Rf for some f e 

2. (a) The moment conditions of Lemma\^ (or Lemma^ are satisfied. 

(b) The solution G{x,t) of the interior problem ^ in C 

( which always exists for t > 0) satisfies the condition 



/dG 
-—(x,t)(j)(x)dx = 
dt ^ ' ^ 

B 



for any eigenfunction (j){x) of the Dirichlet Laplacian in B. 

3. (a) The moment conditions of Lemma\^ (or Lemma^ are satisfied. 

(b) Let — be an eigenvalue of the Dirichlet Laplacian in B and ux be 
the corresponding eigenfunction solution \20\f . Then the following 
orthogonality condition is satisfied: 

j g{x,t)d^ux{x,t)t'''^dxdt = 0. (28) 

5x[0,2] 

Here dy is the exterior normal derivative at the boundary of C. 

4. (a) The moment conditions of Lemma\^ (or Lemma^ are satisfied. 

(b) Letg{x,X) = J g{x,t)jn/2-i{Xt)t^~^dt. Then, for any integer m, 
the m}^ order spherical harmonic term 'gm,{x, A) ofg^x, A) vanishes 
at non-zero zeros of the Bessel function Jm+n/2-iW- 

In fact, if dimension n is odd, one does not need to require the moment 
conditions of Lemma |21 or Lemma |2 since in this case the other range condi- 
tions of Theorem ^1 alone are sufficient. This was first noticed in a different 
setting in 

Theorem 11 Let n > 1 be an odd integer. Then the following four state- 
ments are equivalent: 

1. The function g G C^{S x [0,2]) is representable as Rf for some f G 
C^iB). 
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2. The solution G{x,t) of the interior problem ^ in C (which 
always exists for t > 0) satisfies the condition 

/dG 
-^{x,t)(j){x)dx = 

B 

for any eigenfunction (j){x) of the Dirichlet Laplacian in B. 

3. Let — he an eigenvalue of the Dirichlet Laplacian in B and ux he the 
corresponding eigenfunction solution i2(J\} . Then the following orthog- 
onality condition is satisfied: 

j g{x,t)d^ux{x,t)t'^~^dxdt = 0. (29) 

5x[0,2] 

Here dy is the exterior normal derivative at the boundary of C. 

4- Let^{x,X) = J g{x,t)jn/2-i{Xt)t"'~^dt. Then, for any integer m, the 
mth order spherical harmonic term ^m{x, A) ofg^{x, A) vanishes at non- 
zero zeros of the Bessel function Jm+n/2-i(A). 

As we proved in Lemma El in the case of a ball the polynomial extendibil- 
ity of the i-moments of the data g{x, t) (the moment condition) is equiv- 
alent to the stronger moment condition (Lemma |3)) that requires existence 
of polynomial extensions Qk linked by the additional relations ()25p. 

In fact, if g is in the range g = Rf, f G Co{B), then the polynomials Qk 
obey not only algebraic relations (j25|l . but also the following growth (in k) 
estimates in B: 

mx)\ = \ [ \x~y\''f{y)dy\<2''max\f{y)lxeB. (30) 
J y^B 

B 

It turns out that the moment conditions of Lemma|3]with estimates of the 
above type are not only necessary, but also sufficient for g being in the range 
of the transform R. The following theorem can be regarded as an alternative 
form of Theorem ITUl 

Theorem 12 Let g G C^{S x [0,2]). The following condition is necessary 
and sufficient for the function g being representable as g = Rf for some 
feC^iB): 
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The moments 

oo 





extend from x E S to x E as polynomials Qk{x) satisfying the recurrent 
condition / f^) and the growth estimates 

\Qkix)\<M\xeB, (31) 

for some M > 0. 

The condition of infinite smoothness of functions under consideration is 
not truly necessary. One can prove similar range descriptions in appropriate 
Sobolev spaces, if the functions are supported strictly inside the ball B. 

Theorem 13 Let s > 0. The following four statements are equivalent: 

1. The function g G Hg^^Fiy2(^ ^ ('-'' ^)) '^^ representable as Rf for some 
f^ffcomp^^^^ 

2. (a) The moment conditions of Lemma{^ (or Lemma^ are satisfied. 

(h) The solution G{x,t) of the interior problem ^ in C 

(which always exists for t > 0) satisfies the condition 



/dG 
-^{x,t)(t){x)dx = 



B 

for any eigenfunction (j){x) of the Dirichlet Laplacian in B. 

3. (a) The moment conditions of Lemma{^ (or Lemma^ are satisfied. 

(b) Let — be an eigenvalue of the Dirichlet Laplacian in B and u\ be 
the corresponding eigenfunction solution \20\i . Then the following 
orthogonality condition is satisfied: 

j g{x,t)duux{x,t)t''-^dxdt = 0. (32) 

Sx[0,2] 

Here dy is the exterior normal derivative at the boundary of C . 
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4- (a) The moment conditions of Lemma{^ (or Lemma^ are satisfied. 

(b) Letg{x,\) = f g{x,t)jn/2~i(y^t)t'^~^dt. Then, for any integer m, 
the m}^ order spherical harmonic term 'gm{x, A) ofg{x, A) vanishes 
at non-zero zeros of the Bessel function Jm+n/2-i(A). 

When n is odd, the moment conditions can he dropped. 

Here we used the notation if f (5) for the space of ii<,-functions in the ball 
B with compact support in the open ball. Analogously, iif*^'^^(S' x (0,2)) 
consists of if^-functions on S* x (0, 2) with support in S" x (e, 2 — e) for some 
positive e. 



4 Proof of Theorem [TTH 
4.1 Implication 1^2 

Assume that 1) is satisfied, i.e. g{p,r) = Rf{p,r) for a smooth function / 

supported in B. 

The implication 1 =^ 2(a) is the statement of Lemma |21 

The implication 1 =^ 2{b) is one of the statements of Asgeirsson's theorem 

[T71 1211 EHl El 5 which has already been quoted before. 



4.2 Equivalence 2 <^ 3 

Since conditions 2(a) and 3(a) are the same, we only need to establish the 
equivalence 2(6) <^=> 3(6). This is done in the lemma below. Notice that this 
lemma applies to any bounded domain, not just to a ball. 

Lemma 14 Let D be a bounded domain in with smooth boundary. Ad 
- the Laplacian in D with the Dirichlet boundary conditions, and T > be 
such that all spheres of radius T centered in D do not intersect D. Let also 
geC^idDx[0,T]). 

1. For any eigenfunction (j){x) = (f)k{x) of Ad with the eigenvalue — A^ = 
— A| G a{AD), the following two statements are equivalent: 

(a) 

oo 

dD 
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where 



ux{x,t) = 0(x)j„/2-l(At). 



(b) The solution to the backward initial value boundary value problem 



as t 0+. 

2. If the equivalent conditions (a) and (b) hold for all Dirichlet eigenfunc- 
tions (f) = (pk, then there exists a smooth function f{x) in D, such that 
lim G{x,t) = f{x), and lim Gt{x,t) = 0, where the limits are under- 
stood as convergence in C'^{D) (or, equivalently, in the Sobolev space 
H^{D) for arbitrary s). 

The converse statement also holds, i.e. this behavior of G at t ^ 0+ 
implies (a) and (b) for any eigenfunction 0. 

Proof of the lemma. 

1. Proof of equivalence of conditions 1(a) and 1(b). 

First of all, we need to be sure that a solution G{x,t) of the problem 
exists and is unique and regular for t > 0. This is immediate, due 
to the hyperbolic nature of this problem (at least, until one approaches the 
singularity at t = 0). One can also show this as follows. Applying the 
Weyl transform with respect to time to the functions G{x,t) and g{x,t) in 
one arrives (as we have discussed already) to a similar problem for 
the wave equation, where the corresponding theorems are available in PDE 
textbooks (e.g., j2Hl Section 7.2, Theorem 6]). Then, applying the inverse 
Weyl transform, one obtains the needed solution of (j33|) ^. In fact, a more 
elaborate consideration of this kind can be found further on in this proof. 

^See, e.g., |ST| for usage of such transformation techniques for various PDE problems. 
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We will now prove the implication (a) — (b). We choose a small e > 
and start with a straightforward equality 

T T 

j j gd^uxe-^dSdt = j j„/2_i(At)r-Mt j gd^^ds. 

€ dD e dD 

Using Stokes' formula, one rewrites the inner integral as 



gd^(l)dS = J d^g(j)dS + j {GA(f) - AG0)dx. (34) 

dD dD D 

The first integral on the right is zero, since vanishes on the boundary. Then 
in the second integral, we use the eigenfunction property for (p to get 



gd^(j)dS =- j {AG + X'G)(j)dx. 
an D 

Since G satisfies Darboux equation, we can replace AG by the Bessel operator 
B := -B(„_2)/2 = ^ + This leads to the following form of the last 

expression: 



- J {Gtt + in- l)t-^Gt + X''G)(f)dx. 



D 

Substituting this into the right hand side of (jHll) and changing the order of 
integration, one arrives to 

T T 

gd.uxdSdt =- j <f){x)dx j {Gu + ^^Gt + \^G)jn/2-i{\t)e-'dt. 
e an D e 

(35) 

Let us introduce a temporary notation 

h{t) := / G{x,t)(f){x)dx for t > 0. 




D 



Integrating by parts with respect to t in the inner integral in (jHSj) . we can 



T 



rewrite the resulting expression for j j gd^Uxdsdi as follows: 



dD 




gd^uxdsdt = e"" ^ {ht{e)jn/2-i{et) - h{e){jn/2-i{>^t))t\t=e) ■ (36) 



dD 
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We now need to investigate possible behavior of h{t) and its derivative 
when t — s> 0. In order to do so, let us derive from the Darboux equation for G 
a differential equation for h{t). Applying the Bessel operator d"^ /dt^ + {n — 
l)t~^d/dt to the identity defining the function h{t), one obtains 

'^'4) - - f.OMS. (37) 



df^ t dt J ydt"^ t dt 

D dD 

We used here the Darboux equation for G, integration by parts, the fact that 
(f) is an eigenfunction, and finally the vanishing of (p at dD. 

Let us introduce a shorthand notation for the last integral in (j?fj) : 

w(t) = — J g{x,t)dy(j){x)dS{x). 

dD 

Thus, we get the final non-homogeneous Bessel ODE for h{t): 

d'^h n — ldh^r,^ , , , , 

Due to the condition that g belongs to C^{dD x [0,T]), we conclude that 
w{t) is smooth and vanishes to the infinite order at t = 0. It is a matter of 
simple consideration to show existence of a particular solution of (fHHjl that 
vanishes to the infinite order at the origin. Thus, the type of the behavior 
at the origin is dictated by the solutions of the homogeneous equation. This 
behavior is well known (e.g., ISlj). It depends on whether Bessel functions 
of the first or the second kind are involved. 

Bessel functions of the first kind are smooth and have zero derivative at 
the origin. The ones of the second type, have singularity at zero. If there 
are no Bessel functions of the second kind involved, then the solution of the 
homogeneous equation is continuous at the origin and has zero derivative 
there. On the other hand, if there is a Bessel function of the second kind 
as a part of h{t), then when t 0+, h behaves as follows (e.g., ^^): when 
n = 2, then h{t) = \ogt{C + o{t)) and h'{t) = r\C + o{t)) with non- 
zero constants C. In the case when n > 2, the corresponding behavior 
is h{t) = t2-«(C + o{t)) and h'{t) = ^^-"(C + o{t)). We will now show 
that this type of behavior is impossible, due to (a). Indeed, (a) says that 

/aD gdyUxdsdt — > when e — >• 0. Then, due to 

{ht{e)jn/2-i{et) - h{e){jn,2-i{\t))t\t=e) ^ 0,t ^ 0. 
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On the other hand, if Bessel functions of the second kind were involved, then 
this expression would be C + o(l) with a non-zero constant C, which is a 
contradiction. Thus, we conclude that h{t) is continuous at t = 0, and also 
that h'{0) = 0. The latter statement is exactly the claim of (6). 

Remark 15 In fact, we have proven more than we claimed in (b). Indeed, 
we showed not only that ht(t) = fj^dtG{x,t)(j){x)dx — ^ as t — > 0+, but also 
that h{t) = Jj^G{x,t)(j){x)dx is continuous att = 0. 

The converse implication (6) —>■ (a) is even simpler. Condition (b) means 
that ht{e) — 0, e ^ + . Therefore, h{e) has no singularity at e = and is 
continuous there. Then the right hand side in (jH^ tends to zero as e — >• 
and therefore the left hand side does as well. This means that (a) holds. 

2. Proof of statement 2: Regularity of G at t = 0. In this part of 
the proof, we will use the transformation technique already briefly mentioned 
above, which allows one to toggle between the solutions of the wave equation 
and Darboux equation. Ideologically, what we are about to do, is using the 
Weyl transform. This can be done, and has been done by the authors. How- 
ever, it seemed to the authors, that using only Fourier and Fourier -Bessel 
transforms makes the proof less technical and more transparent. An alter- 
native version of the proof, which uses Weyl transform explicitly is provided 
in Section ini 

First of all, it is well known (e.g., [211 Ch.6.13]) that existence of the 
limit when t — > of G{x,t) and the equality limt^oGt = (even in weaker 
topologies than C°°) mean that G{x,t) can be extended to an even with 
respect to t solution of the Darboux equation. Due to the zero conditions at 
t = T, this even solution will be supported in D x [— T, T]. 

Thus, our task, instead of studying the limits of G and Gt when t —>■ 
(which we would need to do if using Weyl transform), will be to investigate 
existence and regularity (as a function of t with values in H'^{D)) of an even 
with respect to time solution G. 

Suppose we do have such a solution G{x,t) supported in D x [— T, T]. 
Let us then take its Fourier-Bessel transform jFp (with p = {n — 2)/2) with 
respect to time (fT^ . to get a function G{x, A). According to the Lemma [TJ 
this function, as an if*(Z})-valued function of A, would be even with respect 
to A and would satisfy the Paley- Wiener estimate (fTHj) with a = T. Besides, 
the Darboux equation and the boundary conditions would also imply that 
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the following equation and boundary conditions are satisfied: 

{X^ + A,)G{x,\) = OmD ^^^^ 
G{x, A) = g{x, A) for x E S. 

Notice, that ^(x, A) is even with respect to A and of the appropriate Pale- 
Wiener class as a if* (S*) -valued function for any s, due to our conditions on 
smoothness and support of g{x,t). 

So, our problem is now equivalently reformulated as showing existence 
of an even and entire with respect to A solution of (jH^ of the appropriate 
Paley- Wiener class. 

It is clear that (j39j) might not have any solution at all when — A^ belongs 
to the spectrum of the Dirichlet Laplacian in D. However, the necessary and 
sufficient condition for solvability of ()39|) for such values of A are well known 
and easy to derive (they represent the Fredholm alternative): 

J g{x, \)d^(t){x)dS = (40) 
s 

for any eigenfunction of corresponding to the eigenvalue — A^. These 
conditions clearly are equivalent to 3(6) and thus satisfied in our case. Hence, 
one hopes to solve for all A and to eventually get the needed solution. 
This is exactly what we will endeavor now. 

First of all, it will be convenient for us to apply the standard trick of 
moving the inhomogeneity in ()39|) from the boundary condition into the 
equation. Let us denote by E any "nice" extension operator of functions 
from S to D, for instance any one that would map Sobolev spaces H^{S) 
to H^^^^'^{D). Existence of such operators is well known (see, e.g., jb4j). 
The Poisson operator of harmonic extension is one of them. Let us denote 
U{x, A) = G{x, A) — Eg{x, A). Then this function solves the problem 

(A2 + A,)f/(x,A) =/(a;,A) in i? 

U{x, X) = for X e S. ^ ' 

Here /(x. A) = — (A^ + Ax)Eg{x, A)) is of the same Paley- Wiener class with 
respect to A, as g{x, A). 

Let us apply to ()41|) the inverse Fourier (rather than Fourier-Bessel) trans- 
form with respect to A (this amounts to applying the Weyl transform to the 
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original Darboux equation) . Then we arrive to the following evolution prob 
lem: 

^Utt{x,t) = AU{x,t) + f{x,t),x eB,te^^ 
Uix,t)Us = 0. 

Here U and / are inverse Fourier transforms from A to t of t/ and /. Function 
/ is even with respect to t, infinitely smooth as H^{D)-vaAued function of t for 
any s, and is supported (due to the Paley- Wiener theorem) in D x [— T, T]. 
Our goal now boils down to proving existence of an even with respect to 
time solution U{x,t) that is smooth as H^{D)-va\ued function of t and is 
supported in D x [— T, T]. If this is done, then taking Fourier transform 
with respect to time first and the inverse Fourier-Bessel transform next, we 
will arrive to the solution G{x, t) we need, which will finish the proof of the 
Lemma. 

Let us consider the following problem: 

' Utt{,x, t) = AU{x, t) + f{x, t),x e B,t> -T 
U{x,t)Us = (43) 
U{x,-T)=Utix,~T)=0. 

According to the standard existence and uniqueness theorems for the wave 
equation (e.g., j2Hl Section 7.2, Theorem 6]), there exists unique (and smooth 
as H^{D)-valued function of t) solution of this problem. Due to the type of 
the initial and boundary conditions we imposed, one can extend the solution 
to all times by assuming that it is zero for t < —T. 

It only remains to prove that U{x, t) vanishes for t > T and is even 
with respect to time. To do so, let us consider a complete orthonormal set 
{4>k{x)} of eigenfunctions of the Dirichlet Laplacian in D and denote by — 
the corresponding eigenvalues. Let us also expand the functions U and / 
into this basis: 

U{x,t) = J2ukit)(j)kix) 

It will be sufficient for our purpose to show that all functions Uk{t) vanish 
for t > T. 

Let us notice that the following initial value problem is satisfied by Ukit): 

'ul + \luu{t) = h{t) 
Ukit) = u'l^it) = for t < -2. 
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Taking Fourier transform (in distribution sense) in (j45|) . we get 

{Xl-X')M^)=h{X). (46) 

Here /^(A) is even and from the Paley- Wiener class corresponding to smooth 
functions with support in [— T, T] . Consider the function 

As it was mentioned above in this proof, conditions 3(6) guarantee that 
/fc(A) vanishes at the points ±Xk- Thus, Vk is entire, even, and by a simple 
estimate, belongs to the same Paley- Wiener class as fk- Thus, it is Fourier 
transform of a smooth even function Vk{t) supported in [— T, T]. Consider the 
difference Wkit) = Ukit) — Vkif). It satisfies then the homogeneous equation 
w'l + X\wk = and zero initial conditions at t = — T. Thus, it is identically 
zero. Hence, Uk = Vk is even and supported in [— T, T] for any k, and thus 
u{x,t) is also even and supported in D x [—T,T]. This finishes the proof 
of the existence of a solution G+{x,t) of the Darboux equation inside the 
cylinder C that agrees with the spherical mean data on S" x M and which 
is even with respect to time, smooth as an H'^{D)-valued function of t, and 
supported in t e [— T, T]. 

What now remains to prove in the lemma, is the converse statement in 
its part 2. This is, however, trivial. Indeed, the strong convergence of G at 
t ^ we have derived clearly implies the statement lb for any eigenfunction. 

This finishes the proof of Lemma El ^ 

4.3 Equivalence 3 <^ 4 

Since conditions 3(a) and 4(a) are the same, it is sufficient to prove equiv- 
alence of 3(6) and 4(6). Implication 3(6) =^ 4(6) is straightforward. Indeed, 
one can choose in (pUj) instead of ipx one of the eigenfunctions (f)m,i intro- 
duced in fl?I|) . provided A 7^ is a zero of the Bessel function Jm+n/2-i- 
In this case, the integral in (j32|l evaluates to be proportional to gm,i{X) = 

J ^(A, 9)Yj^{9)d9. Thus, vanishing of these expressions for all m, I and A as 

s 

described, is equivalent to the condition 4(6). 

The converse implication 4(6) ^ 3(6) follows analogously, if one takes 
into account the completeness of the system of eigenfunctions (t)m,i- 
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4.4 Implication 2 + 3 + 4^1 

Our goal here is, assuming any of the equivalent assumptions 2,3,4 (or a 
combination of those), to show existence of a function f{x) supported in B 
such that the restriction of its spherical mean Radon transform G{x,t) onto 
the lateral boundary S x [0, 2] of the cylinder C coincides with the function 
9- 

Using the Darboux equation reformulation that we have mentioned be- 
fore, this is equivalent to showing existence in M" x [0, cxo) of a solution 
G{x,t) of the Darboux equation Q such that G{x,0) = f{x), Gt{x,0) = 0, 
and G'|sx[o,2] = 9, for a function / supported in B. Then this function / 
would be a pre- image under the spherical mean transform R of the data g. 

Our strategy consists of solving the following sequence of problems: 

• Showing that the solution G+(x, t) of the interior problem Q is even 
and smooth with respect to t on the whole t-axis. This would, in 
particular, provide us with a candidate f{x) = G{x,0),x G B for the 
pre-image. 

• Using rotational invariance, reducing the problem to single spherical 
harmonic terms of g, G, and /. 

• Showing that each such term G™ of G extends to the whole space 
clS db global solution of Darboux equation. 

• Proving that the value G^{x, 0) is supported inside the ball B and 
coincides with the corresponding harmonic term of /. This will 
show that Rf^ = g'^. 

• Now an immediate continuity argument will show that Rf = g, which 
will finish the proof of this implication, and thus of the whole theorem. 

Let us start realizing this program. 
Interior solution G+. 

Conditions 2(b) and 3(b) mean that the equivalent requirements 1(a) 
and 1(b) of Lemma IT^ are satisfied. Then the second claim of this lemma 
guarantees that the interior solution G+(a;,t) can be continued to an even, 
smooth as H^{B)-va\ued function of t solution of Darboux equation in the 
infinite cylinder B x M.. This resolves the first step of our program. 

Projection to the 0{n)- irreducible representations. 
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Each irreducible sub-representation X"* of the representation of the or- 
thogonal group 0{n) on functions on R" by rotations consists of homogeneous 
harmonic polynomials of a fixed degree m (e.g., |^l98j ). Restrictions of the 
elements of to the unit sphere S are spherical harmonics of degree m. 
We denote, as before, by Yf^{6), I = 1, - ■ ■ , d{m),m = 1, ■ ■ ■ an orthonormal 
basis in X^. The orthogonal projection L'^{S) ^ will be denoted by 

pm. 

{V^h){x) = [ h{y)Z^{y)dS{y), 



where is the zonal spherical harmonic of degree m with the pole x (e.g., 
jnni Chapter 4, Section 2]). Since Bessel operator B and Laplace operator 
A = A^. both commute with the action of 0(n), the projection onto X^ 

G"(x,t) = (P™G) (x,t) 

reduces the Darboux equation. I.e., solves the same Darboux equation 
with the zero data for t = 2 and with the boundary data (7™ = V"^g. Clearly, 
we also have (j™'(x,0) = f"^{x). So, let us assume for now that G = G"^, 
g = g^, and / = . 

Since functions y^*", / = 1, ■ ■ ■ , d{m) form an orthonormal basis of X™, 
we have 

d{m) 

G^{x,t) = E gi{r,t)Yr{e), 

Jim) (48) 

9'^{0,t)= E 9iit)Yr{0). 
I 

As we have already seen, the Fourier-Bessel transform takes the solution 
of Darboux equation to a function G'^{x, A) that satisfies the equation 

A,G'"(x,A) = -\^G'^{x,\) 

in B. Due to this and the special form (j48|) of G™, its Fourier-Bessel trans- 
form can be written as 

d{m,) 

G'-ire, A) = j„/2-i+m(Ar)(Ar)- ^ bi{X)Yr{e), (49) 

1=1 

and correspondingly 

d{m) 

no, A) = jn/2-i+m{x)i^r E (50) 

1=1 
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Now observe that the right hand side of (|49|) is defined for all r, not only 
for r < 1 and therefore defines a smooth extension of G"^(r6', A) for r > 1. 
Thus, we can think of G"^{x,X) as smooth in a; = function defined for 
all {x, A) G M" X M. For x E B, due to Lemma Q this function is of the 
Paley- Wiener class in A, being Fourier-Bessel transform in t of the compactly 
supported smooth function G''^{x,t). However, at this stage we do not know 
much about its behavior with respect to A for x ^ B. To gain this knowledge, 
we need some control over the smoothness and growth of the coefficients bi{X), 
which are defined for all A G C 

Computing the Fourier coefficients with respect to the orthonormal basis 
1^™ of spherical harmonics, we obtain 

biiX){Xry''jn/2-i+m{Xr) = j G'^{r9,X)Yr{9)dS{9). (51) 

6»e5 

The functions bi{X) are clearly analytic at any point Aq 7^ 0. Indeed, for any 
such Aq one can choose r < 1 so that jn/2-i+m{Xr) 7^ for A near Aq. Thus, bi 
is analytic near Aq as the ratio of two analytic functions with non-vanishing 
denominator. 

This argument does not work at A = 0. Moreover, smoothness of bi 
at A = is not guaranteed immediately by (jHH), and requires the moment 
condition. Indeed, let us restrict to the boundary S to get 

6KA)(A)'"j„/2_i+„(A) = j T{0A)Yr{9)dS{9). (52) 

As we already established in Lemma El the moment condition is equivalent 
to the integral in the right hand side in (j52j) vanishing at A = to the order 
at least m. On the other hand, the function X^in/2-i+mi.X) in has zero 
of order m at A = 0. Then, dividing by A™j„/2-i+m(A) in we conclude 
that bi{X) is smooth at A = 0. 

Thus, bi is an entire function. Now we need to estimate its growth at 
infinity (looking for Paley- Wiener estimates). Due to the Paley- Wiener class 
estimates that we have for the expression in the right hand side of (j52p and 
known behavior of Bessel functions, it is a standard exercise to show that their 
ratio bi{X) is of the same Paley- Wiener class as the numerator. Indeed, this 
was treated in J3] . The estimate from below for Bessel functions provided in 
[T^ Lemma 6] and its consequent usage there show that this in fact is true^. 

^An alternative proof can be found in Section IHl 
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Extending 

Now we can apply the inverse Fourier- Bessel transform in A to the ex- 
tended function G^lx, A), x G MJ^ in ()49p . We use the same notation G'™(x, t) 
for the obtained function. This is justified, since by the construction this 
function satisfies the Darboux equation and it coincides with the original in- 
terior "m-irreducible" solution G"^{x,t) in the cylinder i? x M. One can also 
observe that, due to the Paley- Wiener Lemma ^ it is smooth with respect 
to t and supported in 5 x [—2, 2]. 

The size of the support of G""(x, 0). 

Using the relation G"^(x, 0) = f"^{x),x G B and applying the inverse 
Fourier- Bessel transform to (|49|) . one finds 

^ Am) \ 

/-(x) = r{r9) = const J (Ar)™j„/2_i+™(Ar) X"~'dX. 

(53) 

We notice that now we can define an extension F^{x) of f^{x) to the whole 
space M" by applying the formula for r > 1: 

°? Mm) \ 

F™(x) = F'"(r^) := const J (Ar)™j„/2-i+m(Ar) 6KA)V'r(^) A^-^rfA, 

(54) 

which, according to known results (e.g., [HSl Ch. IV, Theorem 3.10]), is just 
the inverse n-dimensional Fourier transform of the following function 

d{m) 

H{x) = H{\e) = X^^hi{\)Yi^{e) (55) 

1=1 

in R" , written in polar coordinates x = X6. Consider for a moment only real 
values of A. One sees immediately that, due to the Paley- Wiener estimates 
on bi, function H{x) is smooth outside the origin and decays with all its 
derivatives faster than any power of |a;| = |A|. If we show smoothness at the 
origin, then H will be proven to belong to the Schwartz class. According to 
standard considerations of Radon transform theory [SJEnilS], this function 
is smooth at the origin if and only if for any k the expression 



U=o 



dX^ '"=°^' 
1=1 
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is the restriction to the unit sphere of a homogeneous polynomial of degree 
k with respect to 6. Due to the form of the last expression, this means that 
^-^^^T^-^(O) = for any k such that either k < m or k — m is odd. The case 
< m is obvious, due to smoothness of bi and the presence of the factor A"*. 
Due to the structure of the functions ^^(A) and Bessel functions, discussed 
already, the condition for k — m odd is automatic. 

Hence, H belongs to the Schwartz space. Then its inverse Fourier trans- 
form, which we previously denoted by F'^{x), is in Schwartz space itself. We 
now need to establish that F"^ is supported inside B. 

By its construction, is the value at t = of a global solution of 
the Darboux equation. Since G"^\sxm.+ = g"^, Asgeirsson Theorem [T7| US] 
implies that RF"^ = 9^- In particular, the integrals of over all spheres 
centered inside B and of radii t > 2 are equal to zero. Indeed, such an 
integral over a sphere centered at x G -B of radius t is equal to G"^{x,t), 
which is known to be zero by construction of C". 

Lemma 2.7 in Ch. 1] claims that if a function decays faster than any 
power and its integrals over all spheres surrounding a convex body B are 
equal to zero, then the function is zero outside of B. In our case we do not 
have all such spheres, but only the ones of radii at least 2 and centered in 
B, where B is the unit ball. However, a simple exercise is to check that the 
proof of the cited lemma still holds and thus is supported in B. This 
means that in fact F™(x) is the zero extension of f^{x) outside the ball B. 

Final step: proving Rf = g 

We already have constructed a function / G C^{B) such that for each 
its component corresponding to an irreducible representation X"^, the 
equality Rf "^ = g"^ holds. Since we have expansions / = J2m /"^ 9 ~ 
J2m9"^ converging in any Sobolev space, the equality Rf = g immediately 
follows by continuity. 

This finishes the proof of implication (2) + (3) + (4) =^ (1) and thus com- 
pletes the proof of Theorem QUI □ 



5 Proof of Theorem [TTI 

We assume now that the dimension n is odd. The claim of Theorem ^2 
is that the statement of Theorem El can be proven without using moment 
conditions of Lemma El or Lemma El In other words, we claim that in odd 
dimensions, the moments conditions (2a)-(4a) of Theorem ITUl follow from the 
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equivalent orthogonality conditions (2b)-(4b). Of course, this effects only the 
sufficiency part of this theorem, since we have proven that in any dimension 
the moment conditions are necessary for g to be in the range. Thus, in odd 
dimensions the moment conditions are redundant. 

In order to prove Theorem ^2 we will follow the part of the proof of 
Theorem ITUl where the moment conditions were not used, and then will finish 
the proof avoiding the moment conditions. 

As we have proven in Section 01 the equivalent conditions (2) - (4) of 
Theorem ^2 (which coincide with conditions (2b) - (4b) of Theorem [rUj) . 
imply the existence of a smooth solution G+(x, t) of Darboux equation in the 
solid cylinder C = i? x M, even with respect to t and such that Gj^{x^ t) = 
for \t\ > 2 and G+I^xr = 9- Here g = g{x,t) is the function introduced in 
Theorem II II We emphasize again that the moment conditions were not used 
in this derivation. 

The key point in the proof of Theorem ^2 is the following auxiliary state- 
ment, which in odd dimensions can be proven without using moment condi- 
tions: 

Proposition 16 The function and all its partial derivatives in x vanish 
on the sphere So = {|x| = 1, t = 0} C M" x {0}. 

We will postpone the proof of this Proposition, and will show now that 
it implies the statement of Theorem ITTl 

5.1 Derivation of Theorem 1111 from Proposition 1161 

Let f{x) := G'+(a;,0). This functions is smooth in B and, according to 
Proposition^! vanishes to the infinite order at S* = dB. Hence, the function 
obtained by zero extension of f{x) to the whole R", is smooth. We will use 
the same notation f[x) for this function: 




Let us now define 



G{x,t) = - [ f{x + ty)dS{y) 



xeW 



n 



yes 
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as the spherical mean transform of this function /. Theorem ^2 will be 
proven, if we show that G = inside C. Indeed, then -R/|sx[o,oo) = 9, and 
thus g is in the range. 

The function G solves Darboux equation in the entire space, and therefore 
the difference 

Z{x,t) = G+{x,t) - G{x,t) 

solves this equation in the cylinder B x [0, oo) (in fact, it extends to an even 
solution in _B X M) . 

Our goal is to prove that the solutions G and G+ agree on the cylinder 
G. At the initial moment t = 0, we have for x G B: 

Z(x,0) = /(x)-/(x) = 0, 
Zt{x,0) = G+,tix,0) - Gtix,0) = 0. 

Then Z vanishes inside the corresponding characteristic cone: 

Z{x,t) = 0, (x,t) e i^i := {0 < t < 1 - 

On the other hand, the interior solution G+{x,t), according to its con- 
struction, vanishes for t > 2. The spherical means G{x,t) also vanish for 
l^^l < 1,^ > 2, since then the sphere {y : \y — x\ = t} does not intersect B, 
and hence also the support of /. Thus, the difference Z = G+ — G vanishes 
in the cylinder {|x| < l,t > 2}. By the same dependence domain argument, 
Z vanishes on the backward characteristic cone K2: 

Z{x,t) = 0, (x,t) G K2 := + 1 < t < 2}. 

Thus, the difference Z of the two solutions vanishes on the union of two 
characteristic cones: 

z{x,t) = 0, {x,t) e K1UK2 

that have the common vertex (0, 1) (see Figure |21). 

Notice that the union Ki U K2 of the two cones contains the segment 
X = 0,0 < t < 2. Each point (0,t) of this segment, except the vertex (0, 1), 
belongs to the interior of K = Ki U K2. Since Z is smooth and Z\k = 0, 
all partial derivatives d^Z{0,t) vanish for all t 7^ 1. By smoothness, this is 
also true for the vertex (0, 1). Thus, function Z vanishes, along with all its 
derivatives, on the entire line a; = 0. We claim that this implies that Z = 
for all {x, t) G G. The following lemma does the job: 
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Figure 2: The cones Ki and K2. 



Lemma 17 Let Z{x,t) be an even and compactly supported in t smooth so- 
lution of Darboux equation in B x M. and xq an interior point in B. If 
{d^Z){xo,t) = for any multi-index a and t G M, then Z{x,t) = for all 
{x,t) e B xR. 

Proof Let us apply Fourier-Bessel transform with respect to t to the function 
Z{x,t). Then Darboux equation transforms to Helmholtz equation and the 
resulting function Z{x, A) is an eigenfunction of Laplace operator: 

AZ{x,X) = -X^Z{x,X). 

By the condition of the lemma, Z{x, A) has at xq a zero of infinite order. Since 
eigenfunctions of the Laplace operator A are known to be real-analytic, we 
conclude that Z{-,X) = 0. Taking inverse Fourier-Bessel transform, we get 
Z = 0. a 

This finishes the proof of Theorem ^2 niodulo the proof of Proposition 
IT^ which we provide in the next sub-section. 

5.2 Proof of Proposition 1161 

The goal of this subsection is to prove that, as the Proposition UHl states, 
at the boundary points (x, 0), the infinite order zero of the boundary data 
g{x,t) at t = implies infinite order zero of G{x,t) with respect to x at 
t = 0, |a;| = 1. 

This proof will be close to the proof of Proposition 7 in 'TI]. Following 
we will reduce the initial-boundary value problem for Darboux equation 



34 



to a problem for the one dimensional wave equation. Again, as in j^I], we 
will use separation of variables in polar coordinates and dimension reduction. 
However, we will do this in a somewhat different manner. Besides, since we 
are dealing with Darboux equation rather than with the wave equation, an 
additional integral transform will be required with respect to t. 
We will break the proof into several steps. 



5.2.1 Separation of variables. 

First of all, we decompose the solution G into spherical harmonic parts G"^ 
that belong to the irreducible representations X*" of the rotation group 0{n): 

dim) 

G{x,t) =^G"(x,t) =^^^;(r,t)r'"Y,'"(^),x = r^. 

m m 1=0 

It is clear that it is sufficient to prove the claim of the Proposition for each 
of these components 

d(m) 

G"™(x, t) = J2 ^K^, t)r"l^'™(^), X = re. 

1=0 

Indeed, we know that G belongs to any Sobolev space H^{B) and hence the 
convergence of the spherical harmonics expansion and Sobolev embedding 
theorems will deduce the claim for G from those for each G^. 

d(m) 

So, we will assume from now on that G = = ^ \E'/(r, t)r'"Y'^™(0). 

1=0 

Correspondingly, we also expand the data g and assume that 

d{m) 
1=0 

The Darboux equation for G{x,t) implies that the coefficients \l/(r, t) = 
\l//(r, t) satisfy the following PDE: 

+ ^^dt^ = d',^ + !i^l±^5,vi/. (56) 

t r 
In order to prove the Proposition, it suffices to prove that 



drP 



:i,o) = o,p = o,i. 
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So, we will concentrate now on studying the solution \1/ of (fHHj) . By its 
construction, \l/(r, t) satisfies the following conditions: 

a) \E'(r, t) is smooth end even with respect to t G M and r G [—1, 1], 

b) ^(r,t) = for \t\ > 2, 

c) ^1/(1,^) vanishes at t = to the infinite order. 

5.2.2 Reduction to the one dimensional wave equation. 

Let us set p = "~^+^™ and apply the inverse Poisson transform (see 0) 
with respect to the variable r: 

(o \ (n-l+2m)/2 

This will reduce the right hand side of equation (j56j) to the second r-derivative. 
We also apply Weyl transform W(„_2)/2 ffTUI) in the variable t. As the result, 
we obtain the function 

U{r,t) = (V-2+2H/2,r>^(-2)/2,^) (r,t), 

which, due to the intertwining properties (O, solves the one dimensional 
wave equation 

Utt - Urr = 0. (57) 

We can observe now that this new function U has the following properties: 

1) t/(±r,±t) = f/(r,t), 

2) U{r,t) = for \t\ > 2. 

3) U{l,t) = q{t) := (P-^W^) (l,t) and q{t) = 0, \t\ > 2. 

The evenness property 1) follows from the evenness of \1/ and the fact 
that both transforms we applied preserve it. Property 2) follows from the 
analogous property b) of ^ and the properties of the Weyl transform W. 
Now, 3) is just the definition of the boundary values of U for r = ±1 combined 
with 2). 

The unique solution of (j37j) in the domain — l<r<l,t>0 with these 
properties, is the sum of the following two progressing waves: 

U{r,t) = q{t-r + l) + q{t + r + l). (58) 

Indeed, this sum clearly satisfies (jFTf) . If r = ±1, t > 0, then f/(±l, t) = q(t) + 
g(t + 2) = g(t), due to 3). So, property 2) holds. Also, q{t — r + l)+q(t+r + l) 



36 



is obviously even in r. Taking into account that q{t) =0 for \t\ > 2, implies 
that q(t — r + 1) + q(t + r + 1) = for t > 2. Thus, by standard uniqueness 
theorem, the solutions U{r, t) and q(t — r + 1) + q{t + r + 1) coincide. 

The representation (jSHj) has important consequences. The first is given 
in the following 

Lemma 18 For all p G Z+, one has 

dPU{±l,0) = 0. 

Proof Let us show first that the evenness of U with respect to t implies 
infinite order zero at of the boundary value q. Indeed, observe that for any 
odd number p = 2s — 1, one has 

g(2s-i)(_^ + 1) + g(2-i)(^ + 1) = {d^^-^U){r, 0) = 0. (59) 

Substituting r = 1 along with using 3) leads to 

g(2^-i)(0) = g(2--i)(0) + g(2^-i)(2) = 0. 

On the other hand, differentiation of with respect to r implies 

-q'~^'H-r + l) + q'~^'\r + l)=0. 

Again, substituting r = 1 yields 

-g(2-)(0) = -g(2-)(0) + g('^)(2) = 0. 

Thus, all derivatives g^^^(O) vanish. This, in turn, implies that U{r,t) has a 
zero of infinite order (with respect to r) at (1,0): 

dPU{l, 0) = (-l)Pg(P)(0) + q^P\2) = 0. (60) 

Since U is even in r, also d^U{—l,0) = 0. □ 

5.2.3 Proving that vanishes to infinite order at |x| = l,t = 0. 

We can now finish the proof of the Proposition^Jby showing that (9^\I'(1, 0) = 
for any p = 0, 1, ■ ■ ■ . Notice, that it suffices to check this identity only for 
even p. Indeed, \l/(r, t) vanishes at (1,0) with all derivatives with respect to 
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t. Then the equation (j56|) imphes that all iterates of the Darboux operator 
acting in the variable r vanish: 

+ ((n - 1 + 2m)/2)drf^{l, 0) = 0. 

Now vanishing of all even order derivatives (9^^\E'(1, 0) clearly implies vanish- 
ing of the derivatives of odd order as well. 

It will be convenient to use the following simple relation: 

Lemma 19 For any smooth even function v{t) the following relation holds: 

\d{t^)J ^ ' \2tdtJ ^ ^ (2j)! ^ ' 

Proof This equality follows from the Taylor formula. □ 

Let us write now the relation ^ = {VW^^)U, using (jH} and (fTn|l . explic- 
itly: 

/ (9 \ (""-^^/^ 

[W)J U{fir,t){l-^i'Y-~'^/'dfi. (61) 

-1 

According to the property c) of the function \E', its all t-derivatives at (1, 0) 
are equal to zero. The odd order derivatives vanish due to the evenness of 
so only the even order derivatives carry interesting information for us: 

9,'^'^(l,0) = 0,p = 0,l,---. 

Let us translate, using (jHH), these equalities into the language of function U. 
Differentiating 2j times with respect to t at the point r = 1, t = under the 
sign of the integral in (jUT|) leads to the expression df-' dj_2~^^^'^U {jj,, 0) . Since 
U{fi,t) is even with respect to t, Lemma IT^ gives 

Taking into account that 2j + n — 1 is an even number, the wave equation 
(or the progressing wave expansion (jiHI)) yields 

9'^+"-^f/(r, 0) = a^^+"-^f/(r, 0). (62) 



38 



Thus, 

1 

^^^■^(1, 0) = const J 9y+"-^f/(/i, 0)(1 - /i2)("-3)/2t;^ = q, (63) 
-1 

for all j = 0, 1, ■ ■ ■ . 

We can now reformulate the identities ()63p as follows: 

Lemma 20 The function B{^) = (9^"^^''^^t/)(/i, 0) is orthogonal in L'^[—l, 1] 
to all polynomials P of degree degP < (n — 3)/2 of the same parity as the 
natural number {n — 3)/2 (i.e. P is even if {n — 3)/ 2 is even, and odd if 
{n — 3)/2 is odd). 

Proof Integration by parts 2j + (n — 3)/2 times in (j63j) and vanishing of 
derivatives of U at (1, 0) lead to 

1 

J B{fi)L^BWfi = 0, (64) 



-1 



where L^ifJ') = const9^(/i^ — 1)"* are Legendre polynomials. Since the deriva- 
tives Lm\j = 0, 1, ■ ■ ■ clearly span the space of polynomials P of degree 
deg P < m that have the same parity as m, this proves the Lemma. □ 

Let us now prove that r-derivatives of vanish at (1, 0). We differentiate 
2k times the identity (jHT|) in the variable r. The wave equation (jFfjl and 
Lemma ^1 imply 

{dfdlr^^^^U){iJ,0) = const(9f +"-i[/)(/i,0). 

Then (jHH) leads to 

^(1, 0) = const / (9f 0)/^''(l - fi^Y'^-^^^^dfx 

-1 

= const / (9f +"-if/)(/x, 0)iJ^^{l - /i2)("-3)/2ci;x. 
-1 

Integration by parts 2k + {n — 3)/2 times gives 

1 1 
^(1,0) = const y"(9^"+^)/2f/)(/i, 0)P(/i)rf/i = const j B{yi)P{y)dii, 

-1 -1 
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where 

is the polynomial of degree (n — 3)/2 and of the same parity as {n — 3)/2. 
Now Lemma EOl claims that all such integrals are equal to zero. Thus, 
9^\1/(±1, 0) = 0, which finishes the proof of Proposition IT^ □ 



6 Proof of Theorem [T2] 

The necessity was already established in (plj). To prove sufficiency, let us 
check that the conditions 3(a) and 3(b) of Theorem hold, which will 
imply that g is in the range of the transform R. 

First of all, the moment condition 3(a) is obviously weaker than our 
condition. 

To check 3(b), observe that the Fourier- Bessel transform 



g{x,X) = J g{x,t)jn/2-i{Xt)r ^dt, 



which is an entire function of A, expands, according to (jl4j) . into the power 
series with respect to the spectral parameter A: 

P oo 

g{x,X) = J g{x,t)jr,/2-i{Xt)r-'dt = J2CkX"'M, 



k[x). 

k=0 



If we replace here the functions Mk{x) by their extensions Qk{x) in the unit 
ball B, then, due to estimates pip for Qk, the extended series uniformly 
converges in B to a function 

oo 


This function ipx is an eigenfunction of Laplace operator: 

oo oo 

A^A = Y,CkX''^kQk-i = ^C,+ic,+iA2(^+i)Q, = -AVa, 

fc=0 s=0 
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due to the relation 

Cg+iCg+i = —CgCs- 

This relation between the coefficients can be easily verified using their explicit 
values 

Therefore, the function ^\{x) = ^{x, A) extends to the ball B as an eigenfunc- 
tion of Laplace operator with the eigenvalue — A^. Let us now apply Stokes 
formula to ipx and a Dirichlet eigenfunction = 0a- Taking into account 
that gx = ipx on S, we obtain 

J gxd.<PxdS = j {i,xMx - ct>xMx)dV = 0. 

5 B 

This provides the orthogonality condition 3(b) (formula ((221) )• Thus, accord- 
ing to Theorem 10, function g belongs to the range of transform R. □ 



7 Proof of Theorem 131 

First of all, the necessity part of the proof of Theorem El clearly survives 
in the Sobolev case. Thus, we only need to establish that a function g e 
H^+^-i)/2^^ ^ (O'^)) satisfying the range conditions, does belong to the 
range of Rs on iff°"^^(i?). Here the following stability estimate is essential: 

Proposition 21 For any e > Q, there exists a constant such that for any 
f G H^{Bi_^) the following estimate holds: 

C~^\\Rsf\\H0^^^_^^^^{Sxie,2-e)) < 11/1^0(^1-.) < Ce\\Rsf\\H°^^^_^^^^iSxie,2-s))- 

(65) 

Here -Bi_e is the ball of radius 1 — e centered at the origin. 

This proposition combines the known injectivity of the operator Rs on com- 
pactly supported functions (see, e.g., j?]) with the elliptic estimates obtained 
recently in [H3] . Such estimates can also be derived from the known ellipticity 
theorem for the pseudo-differential normal operator R^Rs and FIO results 
of jUj. The ellipticity theorem was obtained in |1T] (see also ^^^02], |HH1)- 
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Let us now have a function g{x,t) G -f^s+(„_i)/2('^ ^ (0)2)) satisfying any 
of the range conditions formulated in the theorem. We can choose a positive 
e such that g G -f^°_,_(„_i)/2('^ ^ (2^? 2 — 2e)). Let us extend g to a. function 
gi{x,y) on x M", radial with respect to y, as follows: 

^i(a;,i/) = g{x, \y\). 

The proof of Theorem and lemmas preceding it show that the range con- 
ditions require the following: if one expands gi into a series of spherical 
harmonic terms with respect to x, and then takes the n-dimensional Fourier 
transform with respect to y of each term, the resulting functions have ze- 
ros at certain prescribed locations and of prescribed multiphcities. Let us 
now consider an even, smooth, compactly supported radial approximation 
ipk of the delta-function in M^, such that ^/'fc — > 5 in distributional sense when 
k oo and that the support of ipk shrinks to {0} when k oo. Then the 
convolution gk{x,y) = ipkijj) *y gi{x,y) for a large is a C^, radial with 
respect to y function. If we now denote gk{x,t) = gk{x,y) for \y\ = t, we 
get for large k a smooth function with support in x [e,2 — e] and such 
that gk g in g & ^s+(n-i)/2('^ ^ (^j 2 — e)). Since the Fourier transform of 
the convolution is the product of Fourier transforms, we see that the range 
conditions (being conditions on zeros of the Fourier transform) survive the 
convolution. Thus, functions gk{x,t) satisfy the range conditions. 

According to Theorem ^[ each of the functions gk can be represented as 
gk = Rsfk with fk G C^(-Bi_e). Then the Proposition 1711 shows existence 
of a limit / = \imk-^oo fk in -f^s(-^i-e) such that g = Rsf- This finishes the 
proof. □ 

8 The case of general domains 

Formulation of the range conditions (2b) and (3b) of Theorem [TUl do not use 
explicitly that the set S of centers is a sphere, and hence that the supports 
of functions under consideration are contained in the ball B. Only the range 
conditions (4b) (the ones involving Bessel functions) explicitly use such ro- 
tational invariance. The reader must have noticed that in fact we proved the 
necessity of conditions (2b) and (3b) of Theorem as well as the range 
conditions of Theorem IT^ for arbitrary domain D. The only caveat is that, 
as we have discussed already, the moment conditions (2a) and (3a) should be 
formulated in terms of Lemma E] only, rather than in terms of Lemma |21 In 
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other words, we have proven some necessary range conditions for the spher- 
ical mean transforms with centers on the boundary of an arbitrary smooth 
domain. In order to formulate them, let us introduce a number T such that 
every sphere centered in D and of radius at least T does not intersect the 
interior of D. We now consider the cylinder C = D x [0,T] and formulate 
the conditions 

G(x,T) = Gt(x,T) = 0. (66) 
The following theorem was also proven while we were proving Theorem 

Theorem 22 Let D C M" he a hounded domain with the smooth houndary 
r. Consider the transform 

g{x, t) = Rrf{x, t) = cj-' j f{x + ty)dS{y),x G T, t > 0. 

Ii/I=i 

Then, if f E C^{D) and g{x,t) = Rrf, the following three range condi- 
tions hold. 

1. (a) The moment conditions of Lemma^are satisfied. 

(h) The solution G{x,t) of the interior prohlem / f^) in C 

(which always exists for t > 0) satisfies the condition 

/dG 
-^{x,t)(j){x)dx = 

B 

for any eigenfunction (j){x) of the Dirichlet Laplacian in D. 

2. (a) The moment conditions of Lemma^ are satisfied. 

(h) Let — he an eigenvalue of the Dirichlet Laplacian in D andu\ he 
the corresponding eigenfunction solution ^2(f^. Then the following 
orthogonality condition is satisfied: 

I (67) 

rx[o,T] 

Here is the exterior normal derivative at the lateral houndary 
ofC. 
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3. The moments 

oo 

Mk{x) = j g{x,t)f''+''-^dt 



extend from x & T to x & M."' as polynomials Qk{x) satisfying the 
recurrency condition and the growth estimates 

\Qk{x)\<M\xeD, 

for some constant M > 0. 

Moreover, range conditions (1) and (2) on a function g are equivalent. 

One can ask whether these conditions are sufficient (together with appropri- 
ate smoothness and support conditions on g) for g being in the range of the 
transform i?r- The authors plan to address this topic in another pubhcation. 

One should also note here that in the case of a general domain D, the 
conditions do not take into account geometry of the domain and can be 
specified further. Namely, let us define the following function in D: 

p{x) = sup \x — y\, X E D. (68) 

yedD 

Then clearly can be replaced by 

G{x,t) = for t> p{x). (69) 



9 Proofs of some lemmas 
9.1 Proof of Lemma [U 

As it has been mentioned in text, this is a known result, so we provide a 
quick sketch of the proof here for reader's convenience. 

Let us prove the necessity of the conditions ffist. Evenness is immediate. 
Let us establish Paley- Wiener estimates. Consider the natural extension of 
the function g{t) to a radial function H{y) = g{\y\) on M*^. Due to smooth- 
ness, evenness, and compactness of support of g, we see that if is a smooth 
function on R" with the support in the ball of radius a. The standard Paley- 
Wiener theorem now claims that the ra-dimensional Fourier transform H{^) 
of H is an entire function on C" with Paley- Wiener estimates analogous to 
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(|18|). Since it is known that $(A) is just the restriction of H{C) to the set 
C, = X6, where ^ is a unit vector in M", we get the required estimate (jl8|) . 

Now we prove the sufficiency. Consider function F{x) = ^{\x\) on M". 
Due to conditions on $, function F{x) is smooth everywhere outside the 
origin and decays with all its derivatives faster than any power of So, 
if we can establish smoothness at the origin, this will mean that F belongs 
to the Schwartz class S. Smoothness at the origin, however, immediately 
follows from the radial nature of F and evenness of $. 

Thus, there exists a radial function H on of the Schwartz class, such 
that its Fourier transform is equal to F. If now we define a function g{t) 
such that H{x) = g{\x\), then g is the function we need. We, however, need 
to establish that g has support in [—a, a]. This is equivalent to H having its 
support in the ball |a;| < a. Consider the standard Radon transform of H: 

x-9=s 

Due to the projection-slice theorem [27| ESI EH EOl [72] , the one dimensional 
Fourier transform ^(A, 6) of the Radon data u{s, 6) with respect to the lin- 
ear variable s coincides (up to a non-zero constant factor) with the Fourier 
transform of H evaluated at the point A^, i.e. with F{\9). Due to the 
Paley- Wiener estimates we have for F{\9) and standard ID Paley-Wiener 
theorem, we conclude that u{s,6) vanishes for any 6 and any \s\ > a. Now, 
since H is of the Schwartz class and its Radon transform vanishes for any 
|s| > a, the "hole theorem" HH [701 [721 implies that H{x) = for \x\ > a. 
This concludes the proof of the lemma. □ 

9.2 A Weyl transform proof of Lemma 1141 

We provide here a modification of a part of the proof of part 2 of Lemma [T^ 
It is based on the same Weyl transformation from Darboux equation to the 
wave equation, with the inhomogeneity moved from the boundary conditions 
to the right hand side. So, after the transform G ^ U = WG, where W 
denotes the Weyl transform with respect to t, we, as before, arrive to the 
proving of evenness with respect to time of the solution of the following 



45 



problem: 

Utt - A,U = -dfW{E) in D X [0, T] 

[/(x, t) = for X G dD (70) 

U{x,T) = Utix,T) = for s G 

Here E{x, t) is a smooth function, even with respect to t, and having zero of 
infinite order at t = 0. 

We need to estabhsh that Ut{x, 0) = 0. 

Now the argument starts to differ somewhat from what we had before. 

Let us take any Dirichlet eigenfunction = 0^ in D. It will be convenient 
to rewrite formula (fTUI) using integration by parts and taking into account 
vanishing of the integrand in a neighborhood of oo : 

oo 

(f/(-,t),0) = = const J (9,G(-,s),0)(s2-t2)("-i)/2rfs, 

t 

where (•, •) denotes the inner product in L'^{D). 
Differentiating with respect to t yields 

oo 

t), 0) = -const(n- l)t J{dsG{-,s),^)is^ -t^Y'^-^^/^ds (71) 

t 

We have proven in part 1 that the conditions (a) and (b) imply that 

lim (9tG(-,t),0) =0. 

Therefore, the integral in the right side hand tends to a finite limit as t — * 0+ 
and, taking into account the presence of the factor t in (f7T|l . we get 

{dtU{-,0),<P) = hm{Uti-,t),<P) = 0. 

Thus, the function dfU {x, 0) is orthogonal to arbitrary Dirichlet eigenfunction 
4>{x) and hence 

dtU{x, 0) = 0,x ED. 

Then, due to the wave equation, all derivatives dfU{x, 0) of odd orders p van- 
ish. Hence, iterates of the differential operator d/d{t^) = {2t)~^dt preserve 
smoothness of U{x,t) at t = 0, as a function with values in H^{D). For- 
mula (jlOp implies then that G = yV~^U is continuous and differentiable at 
t = 0. Moreover, ^(^(x, 0) = 0, since dtG{-,0) is orthogonal to all Dirichlet 
eigenfunctions 0. This finishes the proof of the sufficiency part 2 of Lemma 

M 
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9.3 An alternative proof of estimates on coefficients bi 

We provide an alternative growth estimate derivation for the coefficients bi 
(see equation and considerations after it). We would like to establish 
Paley- Wiener estimates for the functions bi{X), A G C. These estimates can 
be derived from the identity (j5ip by averaging over all < r < 1 the point- 
wise estimates that follow from (j5ip . Namely, let us take the absolute value 
in both sides in (j3T|) raised to a fixed power 7 > 1 (to be specified later) and 
integrate with respect to r from to 1. The triangle inequality leads to 

1 1 

\kiX)\^ J |jn/2-i+™(Ar)nArr™dr < J J \G"'ir, X)\^\Yrm'dSie)dr. 

6»g5 

(72) 

Since G"^{x,t) = for t > 2, the Paley- Wiener estimate 
\G"'{re, A)| < Cn{1 + |A|)-VI'™^I 
holds uniformly with respect to r G [0, 1], G S. Then ()72|) gives: 

\bi{X)\^D{X) < constNil + |A|)-^V^I''"^I, (73) 

where 

1 



D{X)= / \jn/2-i+Mr)\^\^rrdr. 







The integral D{X) satisfies the estimate 

|A|D(A) >Do>0. (74) 
Indeed, change of variable r = |A|^^2; in the integral -D(A) yields 

^(A) = \x\-' j |z™j„/2-i+^(^^)rrf^. 


The following estimate of Bessel functions at 00 is well known: 

C 
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Thus, if we take 7 > 2n/{n — 1), the integral converges: 

00 

j \z^in/2-l+m{z)['dz < OO, 


and therefore |A||D(A)| tends, as |A| —>■ 00, to a finite positive constant. This 
gives us the required estimate ()74j) from below. 

Now substituting (fTijl in (f7H|) and raising both sides of the inequality to 
the reciprocal power I/7 yields: 

\bi{X)\ < constN{l + \X\y^e^\''^^\. 

Finally, combining this estimate for bi{X) with Paley- Wiener estimate for 
jn/2-i+m{Xr), we obtain from (jl^ the needed Paley- Wiener estimate: 

IG'^irO, A)| < constNil + |A|)-^e(''+2)l^-^l. (75) 

10 Final remarks 

• The results of Theorems ^1 - ^1 easily rescale by change of variables 
to the ball B of arbitrary radius, which we leave as a simple exercise 
to the reader. 

• It is necessary to note that the range condition (4) of Theorem El is 
an extension of the two-dimensional one in ^]. However, one notices 
a weaker formulation of the moment conditions than in jl5j . 

The condition (3) is a reformulation of (4), but, unlike (4), it is suitable 
for arbitrary domains. 

Condition (2) is not directly the one of jlSli, but they are definitely of 
the same spirit. The authors of [31] work with the wave equation, while 
we do with Darboux. In fact, except the dimension 3, we work with 
different (albeit related) transforms. 

• The condition 1(b) in Lemma IT^ and thus 2(b) in Theorems ITUl and I13| 
as well as condition 2 in Theorem ^2 can be replaced by much weaker 
ones on the behavior of dtG{x, t) or G{x, t) at t = 0. One only needs to 
ensure that the function h{t) constructed in the proof has singularities 
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at t = that are milder than those of Bessel functions of the second 
kind, and thus has no singularities at all. For instance, one can request 
that 



For the same reasons, the above conditions for dtG{x, t) can be replaced 
by analogous conditions for G{x, t), with t^"" replaced by t^""" forn > 2 
and logt for n = 2. 

• It is proven in Theorem ITTl that . similarly to the results of [HTlj, moment 
conditions are not needed in odd dimensions. We suspect that one 
cannot remove the moment conditions in even dimensions, albeit we do 
not have a convincing argument at this point. One notices, however, 
that Huygens' principle played significant role in the proof of Theorem 



One can also notice that, as it is shown in ^3], moment conditions 
alone do not suffice for the complete description of the range. However, 
Theorem shows that a strengthened version of moment conditions 
does suffice. 

• It is instructive to note that we proved in Section 15.11 the following 
statement: Let G+{x,t) be the interior solution of the Darboux equa- 
tion constructed using the orthogonality conditions inside the cylinder 
G = B X M. with zero data at t = 2. If G'+(a;,0) vanishes with all its 
derivatives at the boundary of B, then it extends to the global solu- 
tion of the Darboux equation G, with the initial data supported in B. 
This holds in any (not necessarily odd, like in Section I^TT|) dimension. 
Thus, the job of the moment conditions is to ensure smooth vanishing 
of the interior solution G+ at the boundary, and as we proved, this is 
automatic in odd dimensions. 

• An interesting interpretation of the range conditions (a) and (b) in 
Theorem comes from a model presented in |4 , where a link is es- 
tablished between the spherical mean Radon transform in M'^ and the 
planar Radon transform of distributions supported on a paraboloid in 

We plan to discuss this interpretation in detail elsewhere (see 
also a general approach to such relations in 




D 
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A natural question is: why do we need to restrict the support of the 
function / to the interior of the surface F of the centers? Cannot the 
range of Rr be reasonably described for compactly supported functions 
with supports reaching beyond the surface F? As it was explained 
in this is not to be expected. Briefly, microlocal arguments of 
the type the ones in jUHl IHIl llU3j show that the range would not be 
closed in reasonable spaces (e.g., in Sobolev scale), which is a natural 
precondition for range descriptions of the kind described in this paper. 

One notices that use of microlocal tools as in ^H] was essentially avoided 
in this text (except of the proof of Theorem IT^ . What replaced it, is 
using properties of the solution of Darboux equation instead. Namely, 
existence of the solution G{x, t) in C, and especially its regular behavior 
at t = did the job. Thus, microlocal analysis was replaced by much 
simpler PDE tools (simple properties of the wave equation and Fourier 
transforms) . 

Proposition|2Ilthat concerns stability estimates [HHl, answers a question 
raised in the remark section of [TH] . 

The range condition (2) of Theorem also provides a reconstruction 
procedure: one goes from g{x,t) to the solution G{x,t) of the reverse 
time initial-boundary value problem for Darboux equation in C, and 
then sets f{x) := G{x,0). A similar effect was mentioned in [3T] con- 
cerning their range conditions, where wave equation replaces Darboux. 

Such a procedure is not that abstract. It essentially corresponds to the 
time reversal reconstruction. Similar consideration lead the authors of 
to their reconstruction formulas in odd dimensions (which has been 
extended recently to all dimensions [201123) and the authors of |lUUj 
to what they called "universal reconstruction formula" in 3D. 

One can notice that the proof of our range theorem also involves im- 
plicitly a reconstruction procedure based on eigenfunction expansions, 
assuming the knowledge of the full spectrum and eigenfunctions of the 
Dirichlet Laplacian in B. This procedure (which is somewhat similar 
to the one of [^) would involve at a certain stage division of analytic 
functions with the denominator having zeros. When the data is in the 
range of the spherical mean transform, the range theorem guarantees 
cancellation of zeros. However, such a procedure would be unstable 
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to implement. A better version of an eigenfunction expansion inver- 
sion procedure, which does not involve unstable divisions, has been 
developed recently by L. Kunyansky [SH]. The knowledge of the whole 
spectrum and eigenfunctions is available only in rare cases (besides 
when the domain is a ball), e.g. for crystallographic domains jTHj IT^. 
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